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Abstract 

In  the  Subduction  Area  of  the  eastern  subtropical  North  Atlantic ,  tritium  and  3  He  measurements  were  made  on  sam¬ 
ples  taken  at  over  200  stations  during  the  period  1991-1993.  Combining  these  data  with  salinity,  oxygen  and  geostrophy 
we  can  compute  isopycnal  diffusivities ,  oxygen  consumption  rates ,  and  absolute  velocities.  The  reference  level  velocities 
are  determined  almost  solely  by  the  3H-3He  age  equations  and  to  an  accuracy  of  better  than  Imm/s.  A  classic  absolute 
velocity  spiral  is  seen.  The  depth  variation  of  computed  vertical  velocities  is  consistent  with  vorticity  conservation,  and 
extrapolates  to  within  errors  of  surface  ekman  pumping  and  subduction  rates.  Isopycnal  diffusivities  are  largely  con¬ 
strained  by  the  salinity  equations,  and  are  approximately  1200  m2/s  at  a  depth  of 300m,  and  decrease  gradually  down¬ 
ward.  Diffusivities  apparently  decrease  toward  the  surface ,  likely  an  artifact  of  either  diapycnal  mixing  or  proximity  to 
the  outcrop,  where  unsteadiness  may  be  important. 

The  oxygen  utilization  rates  (OURs)  determined  by  the  oxygen  equations  show  an  exponentially  decreasing  magni¬ 
tude  with  depth.  On  the  shallowest  horizons,  estimated  rates  based  on  a  simple  regression  between  A.O.U.  and  3H-3He 
age  on  isopycnals  agree  well  with  the  more  complete  analysis  done  here,  but  diverge  toward  higher  values  in  the  deeper 
strata.  This  is  an  artifact  of  isopycnal  mixing  combined  with  regionally  concave  A.O.  U.  fields  masquerading  as  consump¬ 
tion  in  the  simplistic  analysis.  Vertical  integration  of  OUR  yields  a  net  water  column  oxygen  demand  of  4. 1±0.8  mol(02V 
m2/y,  which  corresponds  to  an  export  production  of2.5±0.5  mol/m2/y  carbon. 

If  the  3H3He  age  is  a  true  representation  of  the  advective  age  of  fluid  parcels,  then  its  inverse  gradient  for  a  given 
isopycnal  may  be  used  to  compute  the  subduction  velocity  projected  back  to  its  outcrop.  I  present  a  simple  scheme  to  show 
how  to  relate  the  3H-3He  age  to  the  advective  age,  and  after  correcting  for  vortex  stretching,  compute  the  subduction  rate 
as  a  function  of  depth.  The  shallowest  horizon  has  a  subduction  rate  indistinguishable  from  climatological  ekman  pump¬ 
ing  rates,  but  gradually  exceeds  the  projected  rates  with  depth.  The  difference  may  be  attributed  to  the  increasing  impor¬ 
tance  of  buoyancy  forced  subduction  at  higher  latitudes. 


I.  Introduction 


Tritium  (3H)  was  produced  as  a  result  of  atmospheric  nuclear  weapons  testing  in  the  late  1950s  and  early 
1960s,  entering  the  oceans  predominantly  via  precipitation  and  water  vapor  exchange.  It  is  valuable  to  ob¬ 
serve  tritium’s  entry  into  and  movement  through  the  hydrosphere  as  a  dye-like  tracer,  but  its  main  utility  in 
studying  thermocline  ventilation  comes  from  the  fact  that  it  decays  to  a  stable,  inert  tracer:  3He.  Rather  than 
looking  at  the  distributions  of  3H  and  3He  separately,  it  is  useful  to  combine  them  to  compute  a  “tritium- 
helium  age”.  Inasmuch  as  excess  (tritiugenic)  3He  is  lost  to  the  atmosphere  from  the  surface  layer,  and  accu¬ 
mulates  once  a  fluid  parcel  leaves  the  sea  surface,  the  two  tracers  can  be  combined  to  compute  a  3H-3He  age 
that  is  a  measure  of  the  time  elapsed  since  a  water  parcel  was  at  the  sea  surface, 


X  =  X-1log 


1  + 


3He^ 

3H 


With  typical  northern  hemispheric  upper-water  tritium  concentrations  and  present  measurement  technology, 
times  as  short  as  two  or  three  months  or  as  long  as  a  few  decades  can,  in  principle,  be  detected  on  individual 
samples.  Such  time-scales  are  ideal  for  the  study  of  oceanic  subduction,  thermocline  ventilation  and  biolog¬ 
ical  productivity. 

The  purpose  of  this  work  was  to  study  subduction  in  the  eastern  subtropical  North  Atlantic,  and  its  role 
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in  the  ventilation  of  the  main  thermocline.  The  work  involved  obtaining  and  analyzing  tritium  and  helium-3 
samples  from  three  cruises  in  the  Subduction  Area,  and  interpreting  the  results  in  the  context  of 

•  what  are  the  relative  contributions  of  frontal  vs.  distributed  subduction? 

•  can  we  see  the  “smoking  gun”  of  subduction? 

•  what  are  the  rates  of  subduction  as  a  function  of  density  class,  and  are  they  consistent  with  ekman 
pumping  rates? 

•  what  is  the  circulation  of  the  subducted  waters? 

•  what  are  the  oxygen  consumption  rates  in  these  waters? 

•  how  does  the  granularity  of  tracer  distributions  evolve  downstream  (does  it  grow  or  erode)? 

and  to  demonstrate  that  these  observations  can  yield  quantitative  information  about  mixing,  dilution,  ventila¬ 
tion  and  circulation.  We  accomplished  all  of  these  goals  with  considerable  success.  The  highlight  of  this  work 
is  that  we  were  able  to  obtain  estimates  of  subduction  rates  for  the  area,  and  were  able  to  determine  area  av¬ 
eraged  absolute  velocities  to  an  accuracy  of  0.1  cm/s.  In  addition,  we  have  quantitative  estimates  of  vertical 
velocities,  isopycnal  diffusivities,  and  oxygen  consumption  rates. 

2.  Mixing  and  the  Tritium-Helium  Age 

Ideally,  one  would  like  to  interpret  the  tritium-helium  age  in  the  simplest  possible  context:  the  gradient 
in  the  age  is  inversely  proportional  to  the  velocity  on  a  given  surface.  This  might  be  the  case  in  a  simple,  linear 
flow  with  no  mixing,  but  the  tritium-helium  age  must  obey  the  conservation  equation  (Jenkins,  1987) 

pt  =  V(kVx)-u*Vx+1+k(^  +  ^.Vx  (i) 

where  0  is  3H,  £  is  the  sum  of  3H  and  3He,  and  k  is  the  turbulent  diffusivity.  This  equation  resembles  that 
for  “A”,  the  ideal  age  tracer ,  which  obeys 

0  =  V(kVA)-u«  VA+  1  (2) 

except  for  the  unsteady  term  on  the  LHS,  and  the  “pseudo-advective”  term  on  the  far  right.  The  latter  term 
may  be  regarded  as  pseudo-advective  since  it  appears  as  a  scalar  product  with  the  age  gradient  In  younger 
waters  (less  than  a  decade  old)  the  tritium  distributions  are  homogenized,  so  the  resultant  gradients  (and  hence 
this  and  the  unsteady  term)  are  small.  In  older  waters  the  general  sense  of  the  0  and  £  gradients  are  negative 
downstream,  yielding  an  apparent  increase  in  “velocity”  due  to  mixing  non-linearities.  The  magnitudes  of 
these  terms,  in  relation  to  the  others,  are  significant  but  manageable  (see  Jenkins,  1 987  and  below).  The  im¬ 
portant  point  to  note,  however,  is  that  these  terms,  provided  that  k  is  known,  are  all  measurable  quantities, 
and  hence  can  be  accounted  for. 

It  was  expected  that  the  Subduction  Area,  in  the  eastern  subtropical  gyre,  would  be  a  place  were  mixing 
effects  on  the  age  were  small.  This  can  be  seen  with  a  simple,  two-dimensional  gyre  simulation  such  as  the 
one  depicted  in  figure  1.  In  this  simulation,  velocities  were  characteristically  order  1  cm/s  in  the  eastern  part 
of  the  gyre  (consistent  with  the  beta-spiral  results  of  Armi  and  Stommel,  1983)  and  an  isopycnal  diffusivity 
of  1000  m2/s  was  used.  The  entry  of  tritium  and  the  generation  of  3He  on  an  isopycnal  surface  in  an  idealized 
Stommel-type  gyre  were  simulated  as  a  function  of  time  and  "forced”  using  an  extended  version  of  the  Dre- 
isigacker-Roether  (1979)  surface  water  tritium  history,  and  the  tritium-helium  age  computed  from  the  tritium 
and  3He  concentrations.  The  isopycnal  had  an  outcrop  (the  grey  area)  where  some  of  the  streamlines  (indi¬ 
cated  by  the  dashed  lines  in  the  figure)  intersected  the  sea  surface.  Note  that  there  are  two  types  of  streamlines 
in  model.  Open  streamlines ,  which  enter  the  surface  directly  from  the  outcrop  region  represent  fluid  that  is 
subducted  directly  onto  the  isopycnal,  and  hence  may  referred  to  as  the  “directly  ventilated  region”.  On  the 
other  hand,  closed  streamlines  which  do  not  intersect  the  outcrop,  represent  recirculated  waters  within  the 
gyre,  whose  ventilation  must  be  accomplished  by  cross-streamline  mixing.  These  streamlines  form  an  “un¬ 
ventilated  pool  region”  in  the  western  part  of  the  gyre.  The  upper  panel  in  figure  1  shows  the  tritium-helium 
age  distribution  in  the  early  1990s  on  this  isopycnal.  Note  that  the  two  regions  have  a  characteristic  relation¬ 
ship  between  the  isochrons  (contours  of  constant  tritium-helium  age)  and  streamlines:  in  the  directly  venti¬ 
lated  region,  the  isochrons  are  mostly  perpendicular  to  the  streamlines,  whereas  in  the  unventilated  pool,  they 
are  mostly  parallel.  The  former  is  a  signature  of  an  advectively  dominated  system,  while  the  latter  signifies  a 
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mixing  dominated  regime  (Le.,  to  enter  a  streamline,  tritium  must  diffuse  across  streamlines). 

The  lower  panel  in  figure  1  shows  the  100 
influence  of  mixing  on  the  tritium-helium 
age,  by  showing  the  ratio  of  the  tritium-  80 

helium  age  to  the  “ideal  tracer  age”  de¬ 
fined  by  equation  2.  One  can  readily  see  60 

that  in  the  directly  ventilated  region  in  the 
eastern  portion  of  the  gyre,  Le.,  the  Sub-  40 

duction  Area,  the  ratio  is  very  nearly  one, 
indicating  that  the  tritium-helium  age  is  20 

behaving  much  like  an  ideal  age  tracer. 

This  is  consistent  with  arguments  ad-  0 

vanced  earlier  (Jenkins,  1987).  In  the  un¬ 
ventilated  pool  region  in  the  west, 
however,  there  exist  substantial  differ¬ 
ences  between  the  tritium-helium  and 
ideal  tracer  ages,  a  signature  of  mixing.  100 
The  tritium-3He  age  is  generally  de¬ 
pressed  relative  to  the  ideal  tracer  age.  80 

Note  how  iso-ratio  lines  tend  to  be  paral¬ 
lel  to  streamlines  in  the  unventilated  pool  60 

region.  In  general,  though,  the  tritium- 
3He/ideal  age  ratio  is  within  10%  of  1  in  40 

the  eastern  portion  of  the  gyre.  The  con¬ 
clusion  to  be  drawn  from  this  simple  20 

modeling  exercise,  then,  is  that  the  tri¬ 
tium-helium  age  tracer  will  be  dominated  0 

by  advective  effects,  but  some  caution 

must  be  used  to  account  for  the  influence  Figure  1 :  Simulation  of  tritium  penetration  and  the  tritium-helium  age 
of  recirculated  waters.  In  the  analysis  in  an  idealized  Stommel  gyre.  Isochrons  are  solid  lines  (upper  panel) 

which  follows,  this  influence  is  properly  while  streamlines  are  dashed.  Ratio  of  tritium-3He  age  to  ideal  age 

accounted  for.  in  lower  panel. 

3.  Tritium-Helium  Sampling  and  Measurement 

The  Subduction  Area  is  in  the  eastern  subtropical  North  Atlantic,  an  area  expected  to  resemble  the  di¬ 
rectly  ventilated  portion  of  the  thermocline  in  the  simple  gyre  model  described  above.  Here,  one  would  expect 
the  3H-3He  age  to  be  a  more  accurate  measure  of  the  ventilation  age.  As  part  of  this  multidisciplinary,  multi- 
institutional  Advanced  Research  Initiative  by  the  Office  of  Naval  Research,  we  measured  tritium  and  3He  at 
a  large  number  of  locations  and  depths  in  the  eastern  subtropical  North  Atlantic,  between  approximately  40 
and  18°W,  18  and  38°N  from  May  1991  through  May  1993  (figure  2).  Two  cruises,  in  May,  1991  and  May 
1993  were  focused  on  mesoscale  tracer  variability,  and  appear  in  the  map  as  “L”  shaped  arrays.  The  largest 
scale  coverage  was  accomplished  during  two  cruises  in  October-November,  1992.  We  also  obtained  addition¬ 
al  samples  during  December,  1992  taken  for  us  by  T.  Joyce. 

Samples  were  obtained  from  Niskin  bottles  during  standard  CTD-Rosette  casts,  and  transferred  to  90 
cc  valved,  stainless  steel  cylinders.  Dissolved  gases  were  subsequently  extracted  on  shipboard  using  our 
newly  developed  at-sea  extraction  system  into  aluminosilicate  glass  ampoules  (Lott  et  al,  1997).  On  return 
to  the  shore-based  laboratory,  the  gases  were  purified  and  separated  cryogenically  and  helium  isotopes 
measured  using  a  statically  operated,  dual  collecting,  k/2  magnetic  deflection  mass  spectrometer  (Lott 
and  Jenkins,  1984),  while  neon  concentrations  were  determined  using  an  ancillary  quadrupole  mass 
spectrometer.  The  measurement  standards  used  were  volumetric  aliquots  of  marine  air  at  known  pressure, 
temperature  and  relative  humidity,  and  the  results  were  corrected  for  a  small  dependence  of  the  measured 
helium  isotopic  ratio  on  sample  size.  The  running  air  standard  was  repeatedly  compared  to  external  air 
standards  to  ensure  consistency  and  accuracy  The  isotope  measurements  are  reported  as  an  isotope  ratio 


3 


Figure  2:  Tritium-3He  station  locations  (filled  triangles)  for  the  Subduction  Experiment.  The  “L”  shaped  features  are 
densely  sampled  mesoscale  resolving  patterns.  The  solid  circles  designate  the  position  of  current  meter  and  meteo¬ 
rological  moorings. 


anomaly  of  the  unknown  (x)  relative  to  the  atmospheric  standard  (A)  in  per  mil. 


8(3He)  = 


r(3He/4He)x  ^ 

a3He/4He)A  j 


x  1000 


and  are  accurate  to  better  than  1.5  permil  as  determined  by  a  combination  of  standard  reproducibility  and 
comparison  of  replicate  samples. 


Concentrations  of  helium  and  neon,  calculated  using  ion  current  manometry,  again  relative  to  the 
atmospheric  standard,  were  determined  to  better  than  0.2%.  Helium  saturation  anomalies  are  computed 
using  the  solubility  data  of  Weiss  (1970)  according  to 

AHe  =  ^  -  1  j  x  100% 


and  similarly  for  Ne. 

Tritium  samples  were  also  degassed  on  shipboardand  stored  in  vacuo  in  200  cc  aluminosilicate  glass 
flasks.  These  were  returned  to  the  laboratory  and  stored  for  approximately  one  year  for  determination  of 
tritium  by  the  3He  regrowth  technique  (Clarke  et  al ,  1976;  Jenkins  et  al,  1983).  The  tritium  was  typically 
measured  to  an  accuracy  of  0.02  T.U. 


4.  Results 

Figure  3  is  a  map  of  3H-3He  age  on  two  isopycnal  surfaces  in  the  area.  The  shallower  surface  (o0  =26.4 
kg/m3)  outcrops  in  this  area  during  the  winter,  while  the  deeper  (o0  =26.6  kg/m3)  outcrops  further  north.  The 
shallower  surface  is  clearly  tainted  with  freshly  subducted  water,  as  evidenced  by  the  patch  of  water  less  than 
a  year  old  on  the  eastern  portion  of  the  surface.  The  3H~3He  age  increases  to  the  west  and  the  south,  as  would 
generally  be  expected  for  a  southwestward  flow  in  this  portion  of  the  gyre:  the  3H-3He  age  should  increase 
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downstream.  The  deeper  surface  is  characterized  by  both  a  greater  average  3H-3He  age  and  an  apparently 
more  pronounced  age  gradient.  The  “isochrons”  (contours  of  constant  3H-3He  age)  have  organized  them¬ 
selves  in  a  more  orderly  northwest-southeast  orientation.  Note  that  there  is  no  evidence  of  young  (less  than 
one  year  old)  water  on  this  surface,  and  that  the  youngest  ages  are  found  on  the  northernmost  and  shallowest 
portion  of  the  isopycnal  in  the  vicinity  of  the  Azores  Front.  This  suggests  that  the  Front  may  play  a  role  in 
ventilation  of  this  isopycnal  in  the  region,  possibly  by  cross-frontal  mixing  and  lateral  induction  processes 
(e.g.,  see  Qiu  and  Huang,  1995). 

If  the  distributions  of  tritium-helium  age, 
salinity,  and  oxygen  are  known  in  three  dimen¬ 
sions  over  a  portion  of  the  subtropical  gyre,  it  is 
possible  to  combine  them  with  the  field  of  mass 
(via  the  gcostrophic  relation)  to  compute  abso¬ 
lute  velocities  and  isopycnal  diffasivitics,  as 
well  as  oxygen  consumption  rates.  This  is  done 
by  first  computing  the  gcostrophic  velocity  field 
relative  to  an  arbitrary  reference  level  (which 
need  not  be  at  rest),  and  constructing  a  series  of 
conservation  equations  for  those  properties  on 
each  isopycnal.  The  ‘'known”  quantities,  which 
include  the  gcostrophic  velocities,  observed  gra¬ 
dients,  time  derivatives,  etc.  are  substituted  into 
the  equations,  and  organized  so  that  the  un¬ 
knowns  (diffasivitics,  consumption  rates,  and 
reference  level  velocities)  are  on  the  left-hand 
side.  One  thus  obtains  an  over-determined  sys¬ 
tem  of  equations  that  can  be  solved  in  a  least- 
squares  sense  for  the  unknown  variables. 

Samples  were  taken  at  over  200  stations 
during  a  series  of  four  cruises  between  May 
1991  and  May  1993.  Measurement  precisions 
for  3H  and  3He  were  0.02  and  0.03  TU  respec- 
lively,  resulting  in  a  H-  He  age  resolution  of 
approximately  0.25  years  for  an  individual  mea¬ 
surement  and  0.08  years  from  systematic  uncer¬ 
tainties.  The  data  were  interpolated  onto 
isopycnal  surfaces  (in  0.1  kg/m3  increments 
from  =  26.4  to  27. 1  kg/m3)  using  a  simple  local  spline  routine  on  a  station-by-station  basis.  Estimates  of 
large  scale  averages,  isopycnal  gradients  and  laplacians  for  properties,  along  with  attendant  statistical  uncer¬ 
tainties,  are  obtained  by  doing  a  biquadratic  fit  of  the  interpolated  data  on  isopycnals  (solving  the  normal 
equations  using  SVD).  Since  in  general  the  station  spacing  is  comparable  to  decorrelation  scale  of  the  tracers 
(100-1 50km),  the  statistical  uncertainties  associated  with  the  parametric  determinations  is  a  reasonable  mea¬ 
sure  of  uncertainties  in  the  quantities  of  interest.  Laplacians  obtained  from  these  fits  are  thus  averaged  over 
the  whole  region  (an  approximately  1 500  km  square),  but  the  gradients,  evaluated  at  the  center  of  the  region, 
are  a  more  local  average  (perhaps  of  order  500  km  square). 

Time  derivatives  for  the  'H-^He  age  are  determined  by  mapping  the  ^H-^He  age  on  the  same  set  of  isopy¬ 
cnals  for  the  Subduction  Experiment  data  and  for  the  TTO/NATS  (1981)  data  separately,  and  interpolating 
to  a  common  location  (25°  N,  28°  W).  The  age  differences  between  the  two  data  sets  are  then  divided  by  the 
time  differences  to  obtain  the  average  time  derivative  in  x  as  a  function  of  ae  .  Statistics  of  the  fits  suggest 
an  uncertainty  in  this  estimate  ranging  from  0.02  to  0.06,  which  should  be  compared  to  the  age  tendency  term 
of  1 .  The  gcostrophic  calculations  were  performed  using  the  Levitus  ( 1 994)  hydrographic  data-base  interpo¬ 
lated  to  isopycnal  surfaces.  The  gcostrophic  velocities  were  computed  at  the  center  of  the  region  for  each 
isopycnal  from  the  thermal  wind  equation,  with  gradients  computed  using  biquartic  fits  to  the  climatological 
data.  Climatological  data  were  used  to  compute  the  stream  function  rather  than  the  individual  station  data  to 
avoid  transient  eddy  variability  in  the  density  field.  Any  arbitrary  reference  level  may  be  chosen,  and  here  I 


Figure  3:  Tritium-3He  age  distributions  on  the  26.4  (upper) 
and  26.6  (lower)  kg/nr*  isopycnal  surfaces  i*v.  depth.. 
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used  the  1  OOOdb  surface.  For  the  eight  isopycnal  surfaces  in  this  study  ( a0  =  26.4  -  27. 1  kg/m3)  one  obtains 
a  set  of  24  equations,  patterned  after: 


Txlib  +  TyVb- 
for  the  3H-3He  age  (x ), 


-fv2x  +  f^  +  5VvtV  =  l-£ 


at 


xu  —  x  v 

xug  7  g 


(3) 


SxUb  +  Syvb-(V2S)Ki  =  -Sxug-Syvg 
for  salinity  (S),  and 


(4) 


Oxub  +  Oy  vb  -  (V20)Kj  -  Jj  =  -Oxug-Oyvg 


(5) 


for  dissolved  oxygen  (O).  Here  u  and  v  are  the  zonal  and  meridional  velocities,  and  the  subscripts  b  and  g 
correspond  to  reference  level  and  geostrophic  components  respectively.  k{  and  J\  are  the  diffusivities  and 
zero  order  oxygen  utilization  rates,  which  are  assumed  to  be  a  function  of  depth,  and  hence  isopycnal.  Thus, 
one  has  24  equations  in  1 8  unknowns,  which  may  be  solved  by  least  squares.  Prior  to  solution,  the  equations 
were  scaled  (row  normalized)  by  dividing  by  measurement  uncertainties.  The  solutions  (particularly  the  ref¬ 
erence  velocities)  were  insensitive  to  the  precise  scaling  factors. 


4*1  Absolute  velocities  and  Diffusivities 


The  solution  of  the  reference  level  velocity  components  is  dominated  by  the  3H-3He  age  (x)  equations, 
because  the  3H-3He  age  distributions,  especially  on  the  shallowest  isopycnals,  are  dominated  by  advective 
effects.  Further,  the  x  equations  prove  to  be  linearly  independent  for  ub  and  vb  since  the  x  gradient  vector  ro¬ 
tates  significantly  relative  to  the  velocity  vector  with  increasing  depth.  Notably,  the  age  equations  are  mutu¬ 
ally  consistent  within  measurement  errors.  I  can  therefore  solve  for  the  reference  level  velocities,  and  hence 
compute  the  absolute  velocities  for  each  of  the  isopycnal  surfaces.  The  statistical  uncertainty  in  the  magnitude 
of  the  reference  level  velocity  component  thus  estimated  is  0.09  cm/s.  This  may  then  be  recombined  with  the 
geostrophic  velocities,  providing  a  measure  of  the  absolute  velocity  to  an  accuracy  of  0.10  to  0.12  cm/s  for 
individual  isopycnal  surfaces.  The  absolute  velocities,  along  with  estimated  uncertainties  are  plotted  in  figure 
4.  The  velocity  spiral  observed  is  consistent  in  direction  with  that  expected  from  a  P-spiral  calculation  ( cf 
Armi  and  Stommel,  1983)  but  opposite  to  that  originally  obtained  using  a  more  simplistic  treatment  of  3H- 
3He  age  distributions  (Jenkins,  1987).  The  velocities  are  roughly  consistent  in  magnitude  with  the  P-spiral 
velocities  (Armi  and  Stommel,  1983)  obtained  nearby  for  comparable  depth  ranges.  These  velocities  also 
agree  with  the  shallow  Bobber  Float  displacements  (Jim  Price,  priv.  comm.)  but  a  detailed  comparison  is 
made  difficult  by  the  fact  that  the  Floats  track  isothermal  surfaces  that  diverge  from  isopycnals  in  this  area. 

The  isopycnal  diffusivities,  shown  in  figure  5  (left  panel),  are  largely  determined  by  the  salinity  equa¬ 
tions.  Uncertainties  in  the  k  estimates  are  obtained  by  error  propagation  calculations  using  uncertainties  in 
the  salinity  gradients  and  laplacians,  and  in  velocities.  Given  the  relatively  small  size  of  the  mixing  terms  in 
the  3H-3He  age  equations  (figure  5,  right  panel),  uncertainties  in  K  play  a  very  small  role  in  the  determination 
of  ub  and  vb.  The  age  laplacian  term  becomes  increasingly  important  with  depth  due  to  the  combined  effects 
of  deceleration  and  the  rotation  of  the  planetary  scale  3H-3He  age  gradient  relative  to  streamlines. 

The  interesting  feature,  not  discussed  in  the  earlier  p-Triangle  analysis  ( Jenkins ,  1987)  is  the  size  of  the 
unsteady  term.  This  term  is  smallest  on  the  shallowest  surfaces  whose  mean  “age”  is  short  compared  to  the 
time  elapsed  since  the  major  tritium  input,  but  is  larger  for  the  deeper  isopycnals  whose  ages  are  comparable 
to  that  time  frame.  That  is,  the  deeper  isopycnals  are  still  undergoing  significant  tritium  and  £  redistribution, 
and  hence  3H-3He  age  evolution.  The  unsteady  term  thus  becomes  significant  on  the  deeper  isopycnals,  and 
must  be  accounted  for. 

The  oxygen  equations  do  not  influence  the  estimated  diffusivities  (“turning  down”  the  weighting  of  the 
oxygen  equations  by  orders  of  magnitude  does  not  change  the  estimated  k).  The  isopycnal  diffusivities  cal¬ 
culated  are  smaller  on  the  shallowest  isopycnal  and  tend  to  increase  weakly  downward  to  approximately  1200 
m2/s  at  300  m  depth;  and  decrease  very  gradually  below.  Overall  these  diffusivity  estimates  are  lower  than 
that  obtained  for  the  Sargasso  Sea  (Jenkins,  1991),  which  is  consistent  with  large  scale  patterns  of  eddy  ki- 
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Kappa-H  (mA2/s)  Mixing  terms 

Figure  5:  Computed  isopycnal  mixing  coefficients  (left  panel)  and  3H~3He  age  mixing  and  tendency  terms  (right 
panel)  for  the  range  of  isopycnals  from  26.4  to  27. 1  kg/m3. 


netic  energy  {op  cit.).  These  values  are  slightly  larger  than,  but  roughly  consistent  overall  with  values  ob¬ 
tained  from  (3-spiral  calculations  (Armi  and  Stommel,  1983;  Jenkins,  1987)  and  inversions  of  climatological 
data  (e.g.,  Olbers  et  al.,  1985;  Hogg,  1987).  The  deeper  decrease  in  isopycnal  diffusivity  is  consistent  with 
intuition  based  on  the  reduction  of  velocities  (and  shear)  with  depth. 

The  decrease  in  isopycnal  diffusivity  toward  the  surface  is  unexpected,  since  eddy  kinetic  energy  is 
roughly  constant  with  depth,  or  tends  to  increase  upward  {e.g.,  Muller  and  Siedler,  1992).  This  decrease  in  (j 
toward  the  surface  suggests  the  existence  of  unaccounted-for  terms  in  the  salinity  equation  (since  it  dominates 
the  diffusivity  determination)  which  increase  upward  in  importance.  The  two  primary  candidates  are  diapy- 
cnal  diffusion  of  salinity  (KjSzz  )  and  a  possible  unsteady  term  (St ).  The  former  might  arise  from  significant 
vertical  curvature  of  salinity  profiles,  coupled  with  an  upward  increase  in  diapycnic  mixing  rates.  This  can  be 
roughly  evaluated  by  computing  the  magnitude  of  diapycnal  mixing  rate  that  would  be  required  to  become 
significant  relative  to  the  advective  flux  divergence  terms  in  equation  4.  Using  a  number  of  CTD  profiles  from 
the  central  region  of  the  Subduction  area,  I  estimate  the  average  vertical  curvature  in  salinity  at  the  depth  of 
the  Oq  =  26.4  kg/m3  surface  to  be  approximately  5xlO"5  PSU/m2.  Computing  the  isopycnal  advective  salinity 
divergence  to  be  approximately  5x1  Cf 8  PSU/s  using  the  velocities  I  have  obtained  coupled  with  the  salinity 
gradients  obtained  from  the  biquadratic  fits,  I  estimate  that  a  diapycnal  diffusivity  of  order  10'3  m2/s  would 
be  required  to  make  the  terms  equal.  This  value  seems  improbably  large  in  light  of  recent  determinations  of 
diapycnal  mixing  rates  within  the  main  thermocline  (Jenkins,  1980;  Ledwell  et  al.,  1992)  even  with  a  fairly 
strong  depth  dependence. 

Unsteadiness  in  the  salinity  balance  equation  is  most  likely  significant  on  the  shallowest  surfaces,  where 
insufficient  time  has  passed  for  eddy  mixing  processes  to  achieve  the  erasure  of  surface  patchiness  subducted 
with  the  water  mass  and  to  reach  a  statistical  equilibrium  with  the  creation  of  variability  by  the  work  of  eddies 
on  large  scale  mean  property  gradients.  The  time  required  to  achieve  this  may  be  crudely  estimated  from  the 
isopycnal  eddy  mixing  coefficient  (estimated  here  to  be  of  order  1000  m2/s)  and  the  tracer  decorrelation 
length  scale  (estimated  from  tracer  variability  to  be  100-150  km,  see  Jenkins,  1987;  Joyce  and  Jenkins,  1993). 
This  time  scale  is  of  order  1  year  or  less.  Thus  one  would  expect  that  on  the  shallowest  isopycnal,  especially 
those  which  outcrop  locally  and  have  ages  of  order  1-2  years  or  less,  the  unsteady  term  in  the  salinity  equation 
may  be  significant,  affecting  the  shallow  estimates  of  isopycnal  diffusivity.  The  age  and  oxygen  equations, 
incidentally,  are  less  likely  to  be  affected  due  to  the  character  of  their  surface  boundary  conditions  (Jenkins, 
1987;  Doney  and  Jenkins,  1988). 

One  can  also  calculate  the  vertical  velocity  for  each  isopycnal  from  the  absolute  velocities  using  the  re¬ 
lation  w  =  u  •  Vd  where  d  is  the  isopycnal  depth  (right  panel  of  figure  6).  In  addition,  one  can  compute  the 
w  relationship  with  depth  (the  dashed  line  in  figure  6,  right  panel)  using  the  vorticity  conservation  equation 
Pv  =  fwz  and  integrating  upward  from  500m  using  the  tracer  derived  v.  The  computed  w’s  agree  with  the 
vorticity-derived  curve  within  uncertainties  (the  error  bars,  and  the  dotted  curves).  They  also  generally  agree 
with  the  vertical  velocities  obtained  from  the  Bobber  Float  deployments  (the  *  points,  courtesy  Jim  Price). 
For  reference,  the  local  subduction  rate  obtained  from  the  3H-3He  age  data  (the  vertical  gradient  of  the  3H- 
3He  age  may  be  argued  to  be  the  projected  subduction  rate  for  that  density  class)  is  plotted  near  0  m  depth  in 
the  figure  (the  open  circle,  see  section  5).  The  tracer-computed  w  curve  appears  to  extrapolate  to  the  3H-3He 
age  estimated  subduction  rate,  which  also  appears  to  be  consistent  with  the  local  wind  stress  curl  computed 
from  the  ECMWF  gridded  data  (K.  Moyer,  priv.  comm.)  but  a  little  smaller  than  climatological  estimates 
{e.g.,  Isemar  and  Hasse,  1987;  Mayer  and  Weisberg,  1993). 

4.2  Oxygen  Consumption  Rates 

The  distribution  of  oxygen  on  density  surfaces  correlates  negatively  with  the  3H-3He  age,  having  high 
values  in  the  north  and  northeast  and  lower  values  to  the  south.  The  downstream  reduction  of  oxygen  results 
largely  from  in  situ  consumption  due  to  biological  oxidation  of  sinking  particulate  matter  and  possibly  dis¬ 
solved  organic  carbon  {e.g.,  see  Jenkins  and  Wallace,  1992).  One  is  tempted  to  use  the  observed  anticorella- 
tion  between  the  3H-3He  age  and  dissolved  oxygen  on  isopycnal  surfaces  to  directly  estimate  oxygen 
utilization  rates  (Jenkins,  1982,  1987).  However,  the  proximity  to  the  strong  meridional  gradient  in  oxygen 
observed  roughly  along  15  -  20°  N  means  that  one  must  account  for  mixing  effects  in  the  oxygen  distribu¬ 
tions.  Such  effects  lead  to  the  possibility  of  biases  in  such  “determinations”  of  oxygen  consumption  rates,  es¬ 
pecially  in  deeper  waters.  Thus  the  approach  here  is  to  account  fully  for  all  processes  affecting  the  dissolved 
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Figure  6:  Meridional  velocity  component  and  velocity  magnitude  vs.  depth  (left  panel)  and  computed  vertical  ve¬ 
locities  (open  circles,  right  panel).  The  curve  is  computed  from  the  vorticity  conservation  equation,  and  the  dotted 
lines  constitute  1  standard  error.  The  *  are  estimates  based  on  isothermal  float  displacements  (courtesy  J.  Price) 

oxygen  distribution. 

Because  they  appear  only  once  in  each  oxygen  equation  (equation  5),  the  apparent  oxygen  utilization 
rates  (Jj)  are  not  over-determined,  but  are  uniquely  solved  for  within  the  system  of  equations.  Uncertainties 
in  7j  can  be  estimated  by  propagation  of  uncertainties  in  the  velocities,  diffusivities,  and  oxygen  gradient  and 
laplacian  terms.  The  results  are  shown  in  figure  7  (left  panel),  along  with  uncertainties.  For  comparison,  the 
estimated  AOUR  based  on  the  slope  of  the  regression  between  AOU  and  3H-3He  age  for  each  isopycnal  sur¬ 
face  is  included  in  the  plot  (the  solid  line). 

The  oxygen  utilization  rates  (OURs)  in  the  shallowest  layers  are  indistinguishable  within  uncertainties 
from  both  the  regression  based  values,  and  earlier  estimates  (op  cit.)  in  the  p-triangle  region.  This  is  consistent 
with  the  relatively  small  diffusive  flux  divergence  (figure  7,  middle  panel).  The  agreement  is  particularly  re¬ 
markable  in  light  of  the  significant  directional  difference  between  the  oxygen  and  3H-3He  age  gradients  (fig¬ 
ure  7,  right  panel).  Perhaps  this  is  a  further  indication  of  the  unimportance  of  mixing  on  these  horizons. 

OUR  decreases  in  a  quasi-exponential  fashion  with  increasing  depth,  as  would  be  expected  from  an  at¬ 
tenuation  of  the  available  carbon  with  depth,  but  it  does  so  more  aggressively  than  the  regression  based  results 
obtained  both  here  and  from  the  earlier  P-triangle  analyses.  This  increasing  “bias”  with  depth  in  the  regression 
technique  arises  from  increasing  prominence  of  the  diffusive  flux  divergence  with  depth  (see  center  panel  of 
figure  7).  Mixing  driven  by  the  latitudinally  convex  oxygen  distribution,  a  feature  imposed  by  the  larger  scale 
oxygen  distributions,  tends  to  masquerade  as  “consumption”:  as  fluid  parcels  move  down  streamlines,  oxy¬ 
gen  is  mixed  away  (largely  to  the  south),  decreasing  oxygen  concentrations  in  a  fashion  which  resembles  in 
situ  consumption.  Since  the  3H-3He  age  is  not  affected  as  much,  this  biases  the  regression  estimates. 

The  role  of  mixing  and  the  large  scale  gradients  in  the  oxygen  distributions  in  influencing  the  “down¬ 
stream”  change  in  oxygen  also  can  be  seen  in  the  angle  between  the  negative  of  the  oxygen  gradient  and  the 
velocity  vector  (figure  7,  right  panel).  On  the  shallowest  isopycnals,  the  large  scale  negative  meridional  gra¬ 
dient  in  oxygen  is  at  about  a  45°  angle  to  the  velocity,  which  is  oriented  west- south  west  ward.  By  500  m,  the 
velocity  has  rotated  more  westward,  while  the  negative  oxygen  gradient  has  rotated  due  southward.  As  water 
flows  southwestward,  in  addition  to  in  situ  consumption,  oxygen  is  mixed  away  due  to  the  convex  nature  of 
the  large  scale  distribution,  increasing  the  apparent  utilization  rate.  Accounting  for  this  process,  as  has  been 
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done  here,  provides  us  with  a  smaller,  and  more  realistic  consumption  rate  on  the  deeper  isopycnals. 

The  OUR  can  be  least  squares  fit  to  an  exponential  relationship  with  depth  (the  dashed  line  in  the  left 
panel  of  figure  7),  with  attendant  uncertainties  in  the  fit  (dotted  lines).  The  vertical  scale  height  of  the  OUR 
decrease  is  165  m.  Using  this  relationship,  one  can  integrate  the  OUR  vertically  between  120  m  and  550  m 
to  obtain  the  net  oxygen  demand  per  square  meter.  This  integral  is  4.1  ±0.8  mol(02)/m2/y.  Using  the  revised 
02:C  ratio  of  Takahashi  et  al.  (1985),  this  implies  an  export  production  of  2.5±0.5  mol(C)/m2/y.  This  rate  is 
sensitive  to  the  choice  of  the  upper  bound  to  the  integration,  but  appears  to  be  marginally  smaller  than  esti¬ 
mated  productivity  for  the  Sargasso  Sea  (2.8  mol(C)/m2/y,  see  Jenkins  and  Wallace,  1992),  which  is  consis¬ 
tent  with  expectations  based  on  the  availability  of  nitrate.  That  is,  the  regions  of  deeper  mode  water  formation 
in  the  northwestern  gyre  are  areas  of  correspondingly  higher  nutrient  supply  and  hence  greater  new  produc¬ 
tion.  However,  it  should  be  reiterated  that  the  estimated  integral  oxygen  demand  (and  hence  productivity) 
based  on  this  calculation  is  clearly  different  from  the  more  incomplete  treatments  attempted  earlier. 

4.3  Consequences  of  the  3-dimensional  and  larger  scale  characteristics  of  distributions 

There  is  a  substantial  difference  in  the  oxygen  utilization  rate  estimates  between  the  calcula¬ 
tions  presented  here  and  the  less  complete  treatment  of  correlating  AOU  with  3H-3He  age  (see  the  *  profile 
in  figure  7).  The  differences  arise  in  part  from  the  influence  of  the  large  scale  mixing  of  properties  associated, 
for  example,  with  the  confluence  of  waters  of  greatly  differing  origins.  They  also  arise  from  the  complex  in¬ 
terplay  among  mixing,  advection,  and  consumption,  which  must  (and  can)  be  fully  accounted  for. 

The  observed  rotation  of  the  age  gradient  with  respect  to  the  streamlines  as  a  function  of  depth  is 
also  evident  in  the  right  hand  panel  of  figure  7.  On  the  shallower  isopycnals,  the  two  are  virtually  parallel,  but 
proceeding  down  the  water  column,  the  age  gradient  rotates  counterclockwise  (southward)  relative  to  the 
streamlines  that  in  turn  rotate  westward  with  depth.  This  rotation  in  3H-3He  age  gradient  of  course  results 
from  the  increasing  dominance  of  the  planetary  scale  distributions  on  the  age  distribution:  the  same  “front” 
in  properties  seen  in  oxygen  and  salinity  between  15  and  20°  N  is  also  seen  in  tritium  and  age:  this  is  the 
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southern  edge  of  the  £  lens  embedded  in  the  subtropical  main  thermocline.  Also,  the  rotation  is  precisely  why 
the  system  of  equations  used  to  solve  for  the  reference  level  velocity  is  non-degenerate.  The  rotation  of  the 
age  gradient  with  depth  results  in  a  smaller  projected  velocity  u  •  Vx  on  the  deeper  horizons,  and  hence  pro¬ 
duces  an  underestimate  of  the  zonal  component  of  the  velocity.  The  consequence  of  ignoring  this  process  is 
an  apparent  reversal  of  the  tracer  velocity  spiral  (compare  figure  5  to  Jenkins,  1987;  figure  14)  compared  to 
what  is  expected  from  (J-dynamics.  Like  oxygen,  the  presence  of  planetary-scale  property  gradients  plays  a 
significant  role  in  governing  3H-3He  age  distributions. 

4.4.  Subduction  Rates 

At  a  single  location,  the  vertical  gradient  in  the  3H-3He  age  might  be  regarded  as  a  measure  of  the  sub¬ 
duction  rate  for  an  isopycnal  projected  back  to  its  outcrop,  once  corrected  for  planetary  vorticity  gradients 
(see  Jenkins,  1987).  However,  this  approach  assumes  that  the  3H-3He  age  profile  represents  the  true  “advec- 
tive  age”  of  the  isopycnal  surfaces.  The  3H-3He  age  may  unfortunately  be  influenced  by  mixing  processes 
which  in  general  conspire  to  lower  the  apparent  age,  thereby  reducing  the  vertical  age  gradient  and  inflating 
the  apparent  subduction  rate.  This  is  particularly  important  because  the  bias  is  integrated  over  the  trajectory 
of  fluid  parcels.  On  the  shallowest  surfaces,  where  the  age  tendency  term  xt  is  smallest  and  the  extrapolated 
distance  is  shortest,  this  is  likely  not  a  severe  problem,  so  such  an  approach  is  undoubtedly  sound.  Going 
downward  into  the  water  column,  however,  we  are  on  less  firm  footing.  Figure  8  (left  panel)  shows  a  plot  of 
the  inverse  age  gradient  (i.<?.,  l/xz ). 


Inv.  Age  Grad,  (m/y)  Corrected  I.A.G.  (m/y)  Subduction  Rate  (m/y) 

Figure  8:  The  inverse  vertical  age  gradient  (left  panel)  for  pairs  of  isopycnals  between  26.4  and  27.1  kg/m3. 

The  inverse  vertical  gradient  in  the  corrected  age  (corrected  using  date  from  figure  9)  is  shown  in  the  center  panel. 

The  right  panel  shows  the  subduction  rate  computed  from  the  inverse  vertical  corrected  age  gradient  and  the  degree 
of  vortex  stretching  expected  for  movement  from  the  climatological  winter  outcrop  position  of  the  isopycnal  to  the 
center  of  the  survey  region. 

I  obtained  an  estimate  of  the  degree  to  which  the  3H-3He  age  is  representative  of  the  advective  age  by 
running  a  one-dimensional  pipe  model  for  a  variety  of  velocities  and  a  fixed  diffusivity  of  1000  m2/s.It  may 
be  argued  that  this  is  an  oversimplification  of  the  actual  situation,  but  I  would  argue  that  it  captures  the  es¬ 
sence  of  the  tracer  behavior  sufficiently  well  to  make  small  corrections.  The  model  was  run  to  two  points  in 
time  (1979  and  1992)  and  the  advective  age  (defined  as  the  distance  along  the  pipe  divided  by  the  velocity) 
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contoured  as  a  function  of  velocity  and  3H-3He  age.  The  results  are  shown  in  figure  9.  Also  shown  is  the  locus 
of  3H-3He  age  vs.  velocity  for  the  subduction  central  point.  What  is  immediately  apparent  is  that  the  Subduc- 
tion  Area  1992  3H-3He  age  data  are  in  fact  a  good  representation  of  the  advective  ages.  The  3H-3He  age  ac¬ 
tually  overestimates  the  advective  age  slightly,  indicative  of  a  small  overshoot  in  the  3H-3He  age  non-linear 
terms.  I  then  apply  a  correction  to  the  3H-3He  age  data  by  using  the  model  results  to  estimate  the  corrected 
inverse  vertical  age  gradient  (shown  in  figure  8,  center  panel).  Finally,  to  compute  the  isopycnal  subduction 
velocities,  I  further  correct  the  inverse  vertical  age  gradient  for  vortex  stretching  associated  with  movement 
from  the  climatological  winter  outcrop  latitude  to  the  current  location.  This  is  shown  in  figure  8  also  (right 
panel),  along  with  estimated  ekman  pumping  rates  using  annual  averaged  data  of  Hellerman  and  Rosenstien 
(1983)  at  the  winter  outcrop  position  of  the  isopycnal. 

Near  the  surface,  the  ekman  pumping  rate  and  the  subduction  velocity  estimates  agree  within  errors,  as 
would  be  expected.  While  the  ekman  pumping  rate  decreases  with  depth  due  to  the  increasingly  more  north¬ 
erly  outcrop  position  (the  wind  stress  curl  decreases  northward),  the  subduction  velocity  remains  essentially 
constant.  The  difference  between  these  two  curves  may  be  regarded  as  a  measure  of  the  buoyancy  forced  sub¬ 
duction,  which  would  be  expected  to  increase  northward,  and  hence  on  deeper  isopycnals,  due  to  the  geogra¬ 
phy  of  the  winter  mixed  layer  depths  (see  the  arguments  contained  in  Jenkins,  1987).  The  magnitude  of  the 
buoyancy  forced  subduction  rate  required,  however,  is  much  less  than  that  estimated  from  the  p-triangle  data, 
being  of  the  same  order  as  the  ekman  pumping  rates. 

Revisiting  the  subduction  velocity  estimates  made  with  the  1979  data  (Jenkins,  1987)  I  must  conclude  that 
the  deeper  3H--"He  age  data  must  substantially  overestimate  the  true  rate.  This  can  be  seen  from  figure  9  (left 
panel),  particularly  in  the  low  velocity,  high  age  portion  of  the  diagram.  Considering  the  large  magnitude  of 
the  possible  corrections,  it  is  unwise  to  attempt  to  make  a  quantitative  statement  regarding  the  p  -triangle  data, 
but  crude  calculations  indicate  that  the  corrections  indicated  bring  the  estimates  to  within  errors  of  the  Sub¬ 
duction  Area  values  computed  here.  Thus  the  large  disparity  between  the  apparent  subduction  rates  and  ek¬ 
man  pumping  reported  earlier  (Jenkins,  1987)  is  mostly  an  artifact  of  the  non-linear  behavior  of  the  3H-3He 
age  during  the  early  part  of  the  tritium  invasion.  In  1992,  however,  the  tritium  and  3H-3He  age  distributions 
have  now  reached  a  more  steady,  simpler  state. 


1979 


1992 


Figure  9:  Con 
times.  The  sc 


:ours  of  advective  age  as  a  function  of  velocity  and  3H-3He  age  computed  from  a  pipe  model  for  two 
-id  line  marks  the  mean  profile  3H-3He  age  vs.  velocity  from  the  subduction  region. 


5.  Summary  and  Conclusions 
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Consideration  of  a  simple  two  dimensional  gyre  model  suggests  that  one  can  divide  the  gyre  into  two 
regimes:  one  advectively  ventilated,  which  has  streamlines  with  outcrop,  and  one  which  is  mixing  dominated, 
where  streamlines  are  all  closed.  Invasion  of  tracer  into  the  indirectly  ventilated  part  of  the  gyre  must  be  ac¬ 
complished  by  cross-streamline  mixing.  Thus  one  expects  very  diffusive,  low  Peclet  number  behavior  of  trac¬ 
ers  within  that  part  of  the  main  thermocline.  In  the  directly  ventilated  portion  of  the  gyre,  the  3H~3He  age 
behaves  more  as  an  advective  age,  especially  in  the  shallower  layers.  One’s  approach  to  studying  tracer  dis¬ 
tributions  must  be  tuned  to  the  area  of  study. 

In  the  eastern  subtropical  North  Atlantic,  an  area  characterized  by  direct  subduction  and  ventilation,  I 
combined  3H-3He  age  and  other  tracer  data  with  geostrophy  to  calculate  isopycnal  diffusivities,  oxygen  con¬ 
sumption  rates  and  absolute  velocities.  The  reference  level  velocities  were  determined  to  an  accuracy  of  better 
than  1  mm/s.  The  depth  structure  of  the  vertical  velocities  are  consistent  with  vorticity  conservation  and  ex¬ 
trapolates  to  local  ekman  pumping  velocities  at  the  surface.  The  3H-3He  age  equations  dominated  the  deter¬ 
mination  of  the  reference  level  velocity,  while  the  isopycnal  diffusivity  was  largely  determined  by  the  salinity 
equations.  The  isopycnal  diffusivity  is  approximately  1200  m2/s  at  300m  depth,  decreasing  weakly  down¬ 
ward.  The  diffusivity  also  decreases  toward  the  surface,  a  result  which  is  not  intuitive,  and  which  may  be  an 
artifact  of  the  increasing  importance  of  unsteady  or  diapycnal  terms  on  the  shallowest  surfaces.  The  oxygen 
equations  add  no  information  except  in  determining  oxygen  consumption  rates.  Oxygen  consumption  rates 
can  be  integrated  vertically  to  obtain  a  net  water  column  oxygen  demand  of  4.  l±0.8mol(O2)/m2/y,  equivalent 
to  an  export  production  of  2.5±0.5  mol/m2/y  carbon. 

Subduction  rate  as  a  function  of  density  horizon  can  be  estimated  from  the  vertical  gradient  in  3H-3He 
age,  provided  that  it  is  a  measure  of  advective  age.  This  assumption  is  valid  on  the  shallowest  surfaces,  but 
requires  a  small  model-dependent  correction  on  the  deeper  horizons.  The  subduction  velocities  thus  deter¬ 
mined  agree  with  ekman  pumping  rates  near  the  surface,  but  do  not  decrease  with  depth  like  the  isopycnic 
projected  ekman  pumping  velocities.  The  difference  signals  the  increasing  importance  of  buoyancy  forced 
subduction  for  deeper  isopycnals.  The  earlier  (3-triangle  estimates  were  clearly  biased  in  the  deeper  waters 
with  the  non-linear  age  effect,  which  was  stronger  at  that  time.  The  magnitude  of  the  bias  is  estimated  to  be 
sufficient  to  bring  the  p  -triangle  results  into  agreement  with  the  estimates  made  here. 

The  hope  is  that  this  paper  demonstrates  the  utility  and  value  of  transient  tracer  measurements  in  studying 
thermocline  ventilation  and  circulation.  My  thesis  is  that  the  3H-3He  age  and  its  complement  £  provide  pow¬ 
erful  constraints  on  mixing  rates  and  velocities  in  the  subtropical  North  Atlantic.  One  of  the  keys  to  this  suc¬ 
cess  is  combining  the  transient  tracer  observations  with  other  sources  of  information.  Another  is  the 
acquisition  of  appropriate  data  sets.  Clearly  time  series  measurements  are  of  considerable  value  in  document¬ 
ing  the  passage  of  these  tracers  into  the  ocean.  As  important  are  fully  three  dimensional  surveys  of  property 
distributions  to  fully  diagnose  all  aspects  of  mixing  and  circulation.  It  is  extraordinarily  difficult  (if  not  im¬ 
possible)  to  extract  this  kind  of  information  from  sparsely  spaced  sections.  Experience  dictates  that  gridded 
surveys  are  necessary  if  we  are  to  obtain  realistic  estimates  of  oxygen  consumption  rates.  The  work  done  here 
suggests  that  this  approach  can  be  successfully  applied  in  other  subtropical  gyres. 
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